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ABSTRACT 

We identify, categorize, and quantify alignment effects among host and satel- 
lite galaxies using a low-redshift {z < 0.23) sample of spectroscopically-confirmed 
galaxies from the Sloan Digital Sky Survey Data Release 6. Consistent with other 
recent findings, we find that satellite galaxies of red, centrally-concentrated (el- 
liptical) host galaxies with radial velocity separation \AV\ < 600km/s pref- 
erentially reside near the projected major axes of their host galaxies. This 
preference is stronger among red, centrally-concentrated satellite galaxies. We 
explore the dependence of this satellite-host alignment on AV and the pro- 
jected radial separation AR, finding that fractional anisotropy increases with 
decreasing AV and AR. Fractional anisotropy among the closest satellites 
{AR < 250 h^^ kpc) is nearly 40% greater than that seen among most distant 
(SOO/i^^kpc < AR < 1000 /j"^ kpc) companions. We also investigate the ef- 
fects of sample selection and measurement errors on the measured alignment 
signals. Among highly concentrated satellite galaxies at small projected separa- 
tion {AR < 300 kpc), we observe a strong radial (hostward) alignment signal 
in isophotal position angles due to isophotal twisting and contamination that is 
not present when using galaxy model position angles. Among objects for which 
both isophotal and galaxy model position angles agree to within 15°, this elon- 
gation signal is significantly weaker. We also investigate the "Holmberg Effect," 
a well-known result wherein nearby (< AOh"^kpc) satellites of large, inclined 
spiral hosts were seen to preferentially reside near the minor axes of their hosts 
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in projection. Due to the flux limit and spatial sampling limitations, a strict test 
of the "Holmberg Effect" is not possible using only SDSS spectroscopic galax- 
ies. By adopting a looser set of cuts than those of Holmberg's original study, we 
recover a comparable preference of faint blue satellites for the host minor axis 
at marginal (~ Sa) significance. After carefully inspecting the methods used 
by various studies, we suggest that sample selection is largely to blame for dis- 
crepant and contradictory results on galaxy alignment effects in the literature. 
We conclude that several types of alignment likely exist among different galaxy 
populations, but that the observed nature and strength of those alignment trends 
depend sensitively not just on selection criteria but also on the method used to 
determine galaxy orientation. 



1. Introduction 



The study of galaxy alignments is an old and contentious topic. The importance of 
understanding galaxy alignment properties has grown over the past several decades, as ever- 
increasing precision in numerical simulations has revealed more about the galaxy formation 
process. The emergence of weak lensing as a precision cosmology tool has further increased 
the demands on our understanding of the intrinsic alignments of physically-related galaxies. 
An increase in the number of conflicting results merits thorough re-examination of not only 
the existing body of observational results, but also the methodologies used to acquire those 
results, which may be responsible for these apparent contradictions. 

The angular distribution and alignment of satellite galaxies with respect to their hosts 
carries information about the dynamical mechanisms involved in the accretion and subse- 
quent evolution of satellite galaxy halos. For example, it is known from observational and 
theoretical studies that the distribution of galaxies and their host halos is filamentary in 
nature. It has been suggested that large r galaxies accrete satellites along these filaments 
f lKnebe et aU bood : Ihibeskind et al.l boosl : IZentner et all boosi ISailin et al.1 boosh . Subse- 
quent evolution involves complex dynamical interactions between these satellites and the 
host gravitational potential, including dynamical friction and tidal stripping. Luminous 
galaxies within subhalos of larger cluster-size halos allow us to trace the large-scale dark 
matter-dominated cluster potential. The signatures of these dynamical processes may be 
imprinted on the alignments between the central and satellite galaxies. 

While a primary motivation for studying the intrinsic alignment of physically-associated 
galaxies is to provide constraints for galaxy formation models, understanding the magnitude 
of these alignments is important for assessing possible contamination of weak lensing stud- 
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ies (ICroft fc Metzlerl l2000l : iLee fc PenI l200ll ). Since weak lensing of background galaxies 
by foreground substructures induces apparent alignments between these lensed galaxies, in- 
trinsic alignment of galaxies can provide a contaminating signal to weak lensing studies. 
Calibrating the degree of this contamination and its impact on precision cosmology with 



galaxy ahgnment (e.g., 


Croft & Metzler 


2000; 


Lee & Pen 


2001; 


Bernstein & Norbere 


2002 


Mandelbaum et al. 


2005). 



Numerous alignment effects have been observed, many of which are only observed in 
specifi c galaxy popula tions. One of the best known alignment effects, t he "Holmb e rg Ef- 



fect" (iHolmbergi Il969l ). is also currently one of the most contentious. iHolmbergI (119691 ) 
noted that satellite galaxies (SGs) preferentially occupy the space near the projected minor 
axis of very isolated, large, and inclined spiral galaxies. After statistically correcting for 
interloping field objects, there was a significant absence of satellites near the host galaxy 
(HQ) projected major axes. This would seem to indicate that satellites preferentially travel 
or survive along polar orbits. T his result has been both rediscovered at marginal signifi- 



Azzaro et al. 



20o3" in recent years. 



Brainerd 


2005; 


Yane: et al. 


2006 



In the cases of refutation, the inverse effect has been observed. Specifically, these studies 
found that satellites have a stro ng measured prefer ence for the major axis of their host galaxy 
as seen in projection. In both lYang et al.l (120061 ) and lAzzaro et al.l (120071 ). this alignment 
exhibits a color dependence, appearing strongest for satellites of red host galaxi e s. Th ese 
red galaxy populations are obviously dissimilar from those targeted by iHolmbergI (119691 ). 



Galaxy clusters afford another opportunity for alignment analyses. For several decades, 
focused observing campaigns and large- area su rveys al i ke hav e furnished evidence for anisotropy 
in the spatial distribution of cluster members. ISastryl (119681 ) . using a sample of 9 Abell clus- 
ters, provided the first evidence that cluster anisotropy might be common. Over the follow- 



mg years, many more 


1974 




Dressier 




1978 





1981 an d many others). Althou gh 



conflicting results have sporadically emerged (iTucker fc Petersonlll988l : IUlmer et al.lll989l ). it 
is generally agreed that there exists a spatial anisotropy in favor of the projected major axis 
of the brightest cluster galaxy (BCG). This effect is so s trong that the BCG major axis is a 
reliable indicator of the distribution of cluster satellites ( lBinggelilll982l ). This is particularly 
useful when sample flux limits may prevent detection of many cluster members. 

In addition to the alignment signal present within individual clusters, al i gnment has been 



observed between neighboring clusters (lBinggelilll982l ; lArgyres et al.lll986l ; iRhee fc Katgert 



19871 : IWestI Il989al jbl. among others). Although lower in amplitude, this effect is observed 
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to 10 Mpc scales and beyond. Certain recent re sults have linked this alignment to the loca l 
orientation of filamentary large-scale structures (IBailin et al.ll2008l : lFaltenbacher et al.ll2009l ). 



Several studies have also investigated the elong ation of satellite galaxies towards their 
hosts. Numerous studies (e.g.. lCroft &: Metzlerll2000l ) have suggested that an intrinsic align- 
ment of this nature w ould mimic a weak lerising signal and would require careful calibration. 



Othe r recent studies (jPereira fc Kuhrul2005l : lAgustsson fc Brainerdll2006l : iFaltenbacher et al. 



20071 ) find significant evidence of ho stward elongation in an examination of galaxies frorn 



the Sloan Digital Sky Survey (SPSS; kdelman-McCarthv et"aDl200fil: lAbazaii 



Stoughton et al. 2002 : York et al. 200ol ) spectroscopic survey. 



uan 



Faltenbacher et al 



et al 



(2008 



2005 



) find 



a simi lar effect between dark matter halos in N-body simulations. By contrast. iBernstein fc Norber^ 
( 120021 ) observed n o such elongation in a study based on the 2 Degree Field Galaxy Redshift 
Survey (2dFGR,S: [c"olless et al.l[200ll ). 



The size and nature of one's data set plays an important role in the detection and study 
of galaxy alignment effects. Recent large-area surveys (e.g., 2dFGRS & SDSS) contain large 
numbers of galaxies, each with a uniform set of data products. These greater numbers 
allow more stringent selection criteria while still maintaining large statistical samples. The 
uniformity afforded by a fixed reduction pipeline sidesteps most systematic problems with 
intercalibration of different data sets. When available, spectroscopic redshifts can be used 
to eliminate interlopers and greatly simplify selection of physically-related galaxies. 

On the other hand, smaller data sets often benefit from by-eye quality control and selec- 
tion which is impractical for very large surveys. Spurious results that might go undetected 
in an automated pipeline can be found and removed. Some smaller data campaigns also ben- 
efit from fainter limiting magnitudes, which has proven important in studies of individual 
clusters. 

The purpose of the present study is to independently investigate the alignment effects 
mentioned above using the spectroscopic ga laxy catalog of the Sloan Digital Sky Survey Data 
Release 6 (lAdelman-McCarthy et al.ll2008l . hereafter SDSS DR6). We also seek to resolve 
any discrepancies between our results and those that have been presented to date. 

Using a subsample of low-redshift {z < 0.23) galaxies from the SDSS DR6 spectroscopic 
catalog, we explore the distribution and alignment properties of satellite galaxies in projec- 
tion. Groups are determined by centering a cylinder of fixed physical size (2400 km/s in 
"height," 1 Mpc in radius; Fig. [3^) on the brightest local galaxy and counting objects 
that fall within the cylinder volume. We impose a variety of different cuts in observed galaxy 
properties and observe the resultant changes in satellite galaxy distribution and alignment. 



The main focus of this study is quality control in both selection criteria and measure- 
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ments of galaxy properties. To do this we simultaneously use multiple independent measures 
of galaxy position angle and shape from SDSS. To determine reliable relative luminosities, we 
incorporate all 5 photometric bands through spectral templates, improving host determina- 
tion among different galaxy types. Lastly, we examine distribution and alignment properties 
separately as functions of host and satellite color and concentration, host-satellite velocity 
separation, and projected radial separation, combining and improving upon results found in 
the literature. 

The content of the present work is organized as follows. In [|2] we describe the data 
products used in the study. In ^ we describe the selection criteria used to pare the initial 
data set and detect groups within it. In §l]we present the properties of our data sample and 
discuss our methods of quality control. We present our findings in ^ Numerous caveats 
and complications follow in ^|6l How these results compare to previous findings follows in 
^ We present our investigation of the Holmberg Effect in ^ Lastly, a brief summary and 
discussion of our findings and their implications can be found in ^ 

In what follows we adopt the following cosmological parameters: Q/^ = 0.7, Qm = 0.3, 
and h = i/o/(100 km/s/Mpc) = 0.7. 



2. Data Set 



The galaxies we use are a subset of the SDSS DR6 (lAdelman-McCarthy et al.l l2008l ) 
spectroscopic galaxy sample. In its entirety, the main spectroscopic catalog contains 790,220 
galaxies to limiting Petrosian magnitude < 17.77. We include all spectroscopically- 
confirmed galaxies with well-established {zConf> 0.35) redshifts of 0.004 < z < 0.23. This 
lower bound safely excludes possible neighbors of our own galaxy while the upper bound 
minimizes the impact of fiber collisions. We further select galaxies with g > and r > 
(modelMag, see below). Lastly we require successful (not placeholder) measurements of 
isophotal axes {isoA & isob)- The resultant sample contains 572,495 galaxies which meet 
these criteria. 

SDSS employs matched de Vaucouleurs and exponential galaxy models to optimally 
extract many photometric properties, including some that we use in this study. A de Vau- 
couleurs i?^/^ profile accurately describes the li ght profiles of many el liptical galaxies and 



spiral bulges. SDSS employs the following form (IStoughton et al.ll2002l ): 



I{R) = Jo exp [-7.67 {R/R,fff'^]. (1) 

This profile is truncated beyond 7Reff, decreasing smoothly to zero at SRejf. The profile is 
softened slightly within R < Reff/50. 



- 6 - 



Exponential light profiles are frequently used to describe dwarf ellipticals and the disks 
of spiral galaxies. The SDSS exponential model, 

I{R) = Ioexp[-l.Q8{R/R,ff)l (2) 

is truncated beyond SR^jj, smoothly decreasing to zero at 4_Reif fjStoughton et al.ll2002l ). 



The SDSS pipeline convolves two-dimensional variants of the above profiles with a 
double- Gaussian approximation of the locally-determined PSF in order to account for the 
effects of seeing. The best fit of the exponential and deVaucouleurs models to each galaxy 
simultaneously determine the axis ratio, position angle (PA, in degrees East from North), 
and effective scale radius {Rejj^ in arcseconds) in each case. The axis ratio and position angle 
are taken to be constant throughout the model galaxy. Although fitted ac ross a range of 



radius, these models tend to best trace the bright, inner regions of galaxies (IStrateva et al. 
2nnil : Istoughton et al.l[2nn^ . 



SDSS provides several different galaxy flux measurements. Unbiased colors, which are 
necessary to determine K-corrections to the highest possible accuracy, must be obtained 
using the same aperture in each wavelength bands. To accomplish this we employ modelMag 
(model magnitude, hereafter simply magnitude) fluxes in the present study. Both matched 
galaxy models (see above) are fitted to each object. The fit of higher l ikelihood in r-band 
is re-apphed to w, q, r, and allowing only the amplitude to vary (ILupton et all 12001 



Stoughton et al.lbood l The resultant measurement, termed modelMag, thus uses the same 
effective aperture in all bands. 

SDSS additionally computes a best-fit composite model flux in each band. This cmodel 
magnitude very accurately captures the total flux of each galaxy by calculating the best-fit 
non-negative linear combination of exponential and de Vaucouleurs profiles separately for 
each of the five photometric bands. The fractional flux contribution of the de Vaucouleurs 
profile, fracDe V, is a useful (seeing-corrected) measure of light profile concentration. 

The Petrosian fiux, also computed for each galaxy, measures the light within a circular 
aperture the radius of which is determined by the slope of the azimuthally-averaged galaxy 
brightness profile. In t he SDSS, this Petrosian radius (rp) is chosen such that 7^p(rp) = 0.2 



(Istoughton et al.ll2002l ). where TZp{r), the Petrosian ratio, is the ratio of sur face brightness 



in ail annulus at radial position r to the mean su rface brightness interior to r (IBlanton et al. 
200ll : lYasuda et allboOll : Istoughton eraDl2002h : 

f'f'rfr' 27rr7(r')/[7r(1.252 „.„ ^ 

The Petrosian fiux is then computed by summing the fiux within 2rp. As implemented by 
SDSS, Tp measures the size of the bright inner regions of galaxies. Its distribution closely 
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matches (without artifacts) those of the exponential and de Vaucouleurs model radii. Using 
Tp we reject extremely close neighbors that are likely heavily blended or misidentified bright 
knots within a larger galaxy. In our study we exclusively use the r-band measurement. 

Galaxy light profile concentration is frequently measured using the inverse concentration 
ratio (C), defined as the quotient of Petrosian 50% and 90% radii \C = R^o/Rm, where Rr,n 
a nd -Rqn are the radii that contain 50% and 90% of Petrosian flux (IStoughton et al.ll2002l : 



Blanton et al 



2001 



(IShimasaku et al. 



2001 



Yasuda et al.l I2OOII)]. C is highly correlated with morphological type 



Strateva et al.ll200ll ) and provides a convenient and effective way to 



statistically separate large galaxies into different morphological classes. The Petrosian radii 
are not, however, corrected for the effects of seeing, reducing the accuracy of concentration 
ratio for smaller objects of sizes nearing that of the PSF. Furthermore, as seeing is depending 
on the time of observation, relying on C alone may introduce time-dependent systematic 
effects. 

In the present work, we measure light profile concentration with SDSS galaxy parameter 
fracDeV. Unlike Petrosian radii R^q and Rgo, the galaxy models are corrected for seeing 
conditions (see above) and should thus provide a straightforward and reliable measure of 
concentration which will depend less on the angular scale of the individual galaxy relative to 
that of the PSF. We illustrate the distributions of color and inverse concentration and how 
they relate to fracDe V in Figure [H 

Although the galaxy fluxes are very well characterized, systematic effects in the model- 
fitting procedure lead to some discreti zation in the determined model scale radii and con- 
sequently into the model axis ratio (lAdelman-McCarthy et al.l |2008| . see §6.11) . There- 
fore, to quantify the gal a xy shape, we instead calculate second-order adaptive moments 
(IBernstein fc JarvisI l2002l : iHirata fc Seljakll2003l ). The distribution of galaxy shapes deter- 
mined in this way closely matches the distribution determined by fitting the de Vaucouleurs 
and exponential models, but the adaptive moment method is free of the model-fitting arti- 
facts. 

Adaptive moments shape measures are easily calculated using data provided by SDSS. 
SDSS DR6 uses a Gaussian weight function w{x,y) matched to the shape and size of the 
galaxy image I{x,y). The first-order moments 

^ /xw(x,y)/(x,y) dx dy 
J w{x, y) I{x, y) dx dy 

yield the galaxy positions. Using this value, the second-order adaptive moments can be 
computed as in 

f{x — xq)"^ w{x, y) I{x, y) dx dy 



J w{x, y) I{x, y) dx dy 



(5) 
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The SDSS DR6 provides the second-order parameters r, e+ and Cx for all galaxies, where 
r = M^x + e+ = {Mxx — ^yy)!^-, a^d Cx = 2Mxy/T. Using these values, we define the 
adaptive moments axis ratio 

Qmom = , (6) 



where e = (e^ + e^)^^^ ( Rydenll2004l ). This shape is not corrected for effects of seeing. For 
objects larger than the PSF (true of the objects in our sample), however, any such correction 
will be negligible. In our analysis we exclusively adopt the value of qmom determined from 
the r-band measurements. 

SDSS provides model-independent measurements of PA and ellipticity by fitting to the 
25 magnitudes / arcsec^ isophote. After first determining the object centroid, the reduction 
pipeline determines the radius of this isophote as a function of angle on the sky. This 
angular profile is then expanded in a Fourier series, the coefficients of which determine the 
isophotal major {isoa) and minor {isob) axes, plus position angle. Although this procedure 
is performed separately in all 5 photometric bands, we use the r-band PA exclusively for 
both isophotal and galaxy model PA. 

Although anticipated in a future Data Release, SDSS does not presently provide error 
estimates for its position angle (PA) measurements. We note, however, that the galaxies in 
the spectroscopic sample are considerably brighter than the SDSS photometric limit. As a 
result, for those objects that meet our shape criteria (see §3.ip . PA measurement errors should 
be negligible. We do, however, observe significant systematic errors in the PA measurements. 
These must therefore be treated with care (see §6.ip . 

From isoA and isob, we compute both the axis ratio and an isophotal "radius." We 
define isophotal axis ratio qiso = isoA/isos (0 < qiso < 1) and isophotal radius isor = 
y/isoA isoB- We use the latter to exclude spurious satellite galaxies (e.g., blended objects or 
bright substructure) and to eliminate possible cases of isophotal contamination (see §6.11) . 

In practice, the 25 mag arcsec"^ isophote tends to trace the outermost observable re- 
gions of galaxies. The measured isophotal PA may, in fact, be very different from the PA 
determined from fitting models to the galaxy light distribution. In our analysis, we compute 
the position angle difference (APA) between isophotal and model PAs which we then use to 
identify cases where one (or both) PA measurements are unr eliable or where the two PA mea- 



surements are genuinely discrepant. Several recent studies f lBrainerdll2005l : lYang et al.ll2006 



Faltenbacher et al.ll2007l ) have reported alignment findings based on isophotal PA alone. We 
believe that these two different PA measurements can have different physical meaning with 
respect to galaxy structure. As a result, both should be considered. We explore the effects 
of using one or the other PA determination in §6.1[ 
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FoUowing lBailin et al.l (120081 ). we K-correct our galaxy colors to z = 0.1. Selection of this 
"neutral" redshift {z = 0.1 is both the mean and median of redshift our sample) minimizes 
K-correction error for the bulk of our galaxies. We denote K-corrected magnitudes in the 



usual way with the redshift in superscript (e.g., " r). To p erform the correcti ons, we employ 
the Low- Resolution Template^ (LRT) software package of Assef et al. ( 2008 ). This package 
also uses these templates to calculate the bolometric luminosity (where Lboi spans 0.2 to 10 
fim), which we use to robustly separate host and satellite galaxies. 

Finally, we use this K-corrected phot ometry to assign ea ch galaxy one of two color types 
(red or blue) using the following criteria ( Bailin et al. 2008 1^: 



red galaxies: 
blue galaxies: 



r) > 0.78 - 0.0325(M^ - 5 log /i + 19) (7) 
r) < 0.78 - 0.0325(M^ - 5 log /i + 19) (8) 



This effectively separates the so-called red sequence and blue cloud (Fig. [2]). These data 
products allow us to isolate and characterize the populations which give rise to the alignment 
trends we observe. 



3. Group Selection and Alignment Geometry 

Before analyzing alignment effects, we must first identify physically-related galaxies. To 
do so we separate galaxies into groups. We define a fixed-size cylindrical "neighbor volume" 
of 2400 km/s in height and 1 h^^Mpc in radius (Fig. [3^) centered on each galaxy, within 
which we search for companions. Each group consists of one host galaxy (HG) and at least 
one satellite galaxy (SG). 

Initially, any galaxy in our data set is a potential HG. We begin the selection process 
by eliminating all galaxies whose neighbor volume intersects either an SDSS spectroscopic 
survey border or the limiting redshift of our sample. For each remaining object, we produce 
a list of companions (all galaxies within the neighbor volume), eliminating any solitary hosts. 
Next, using the bolometric luminosity calculated with the LRT code (||2]), we eliminate all 
galaxies which have a more luminous companion within the neighbor volume. The remaining 
objects comprise our final HG list. 



^The Low- Resolution Templates of ^ Assef et al. ( 2008 ) are a collection of publicly available Fortran rou- 
tines. See http : //www . astronomy . ohio- state . ed.u/"^r j assef /Irt for more information. 

^This division deviates slightly from that presented in Bailin et al. ( 2008h . We use the values from version 
1 of their arXiv preprint ([Bailin et alj|2007| ). The differences between the two divisions are minor and do 
not affect our results. 
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Our group sample consists of all companions (SGs) within the neighbor volume of all 
HG list members. In total, we find 291,435 SGs in the vicinity of 101,850 unique HGs. In 
our implementation, SGs can belong to multiple groups while the HGs cannot. 



3.1. Alignment Angles 

Several distinct alignment trends have been observed to date (see Identifying and 
quantifying these effects requires more than one angular separation measurement. For each 
HG-SG pair, we compute a location angle (6*) and radial alignment angle (0). These quanti- 
ties characteri ze the geometries w'e inves tigate (see Fig. [3]). For simplicity we have adopted 



the symbols of iFaltenbacher et al.l (120071 ). We separately calculate each using isophotal, ex- 
ponential, and de Vaucouleurs PA values. To exploit inherent symmetries, we reduce all 
angles to values between ° and 90 °. 

We define 9 to be the angle between the HG major axis and the line connecting the 
centers of host and satellite (i.e., the host-satellite separation vector). This quantity describes 
the location of a companion galaxy with respect to the major (^ = 0°) or minor [9 = 90°) 
axis of its host (see Fig. [3]). We use this to characterize anisotropy in the satellite galaxy 
distribution (i.e., the Holmberg Effect). 

Calculating 6 accurately requires a robust measurement of the HG PA. To this end, we 
include only those HGs with qiso < 0.9 and qmom < 0.9. This shape requirement eliminates 
78,070 galaxy pairs (27%). We also exclude all objects with APA > 15°. By itself, this 
cut eliminates 90,341 pairs (31%). In combination, the shape and APA cuts eliminate 43% 
of our group sample but ensure that the remaining galaxy pairs (166654 SGs around 58259 
unique HGs) have usefully high HG PA accuracy and thus reliable 6 measurements. This 
procedure maximizes the sample size without risking inclusion of spurious PA measurements. 

Our second angular separation, 0, is the projected angle between satellite galaxy position 
angle and the projected host-satellite separation vector (Fig. [3]). We use </> to expose any 
preferred direction in SG elongation. The orientation extrema are radial (0 = 0°, directly 
toward the host) and tangential (0 = 90°). 

To properly measure 0, we require a robust measurement of SG PA. We thus include 
only those SGs with qiso < 0.9 and qmom < 0.9. Applied to SGs, this shape cut ehminates 
66,159 pairs (23%). We also eliminate SGs that lack good agreement between model and 
isophotal PA (i.e. SG APA > 15°). On its own, this APA requirement eliminates 85,105 
galaxy pairs (29%). Together, these criteria reduce our sample size by 38% (to 180,777 SGs 
around 77,213 unique HGs) but minimize the inclusion of spurious PA measurements. We 
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do not apply any shape cuts to HGs in measuring (j). 

For each object, we calculate angular separations using spherical trigonometry, adopting 
the galaxy centers indicated by the SDSS J2000 coordinates. Position uncertainty is negli- 
gible for our purposes. SDSS determines centroids (the first moment of light di stribution) 



using an adaptively-smoothed, PSF-length-scaled quartic interpolation algorithm (jPier et al. 



20031 ) which is accurate to within a few tens of milliarcseconds. The greatest sources of un- 
certainty in calculating 6 and are due to the PA measurements. We follow the shape and 
APA prescriptions outlined above to minimize the potential impact of ambiguous, uncertain, 
or contaminated PA determinations on our results. 



4. Statistical Sample Properties 

Before the cuts on HG ellipticity and position angles, the group sample consists of 
291,435 satellite galaxies (SGs) within the cylindrical neighbor volume of 101,850 unique 
host galaxies (HGs), indicating a mean SG count of 2.86 per system. Nearly half of systems 
(48,018) are galaxy pairs and 99% of systems contain 20 or fewer SGs. The largest observed 
group contains 169 satellites (Abell 2197). 

In FigureHlwe present spatial, velocity, concentration, and shape distributions of objects 
in our sample of grouped galaxies. Figure shows the number of SGs within annular bins 
of projected radial separation {r'^dN/dr, where r is the HG-SG projected radial separation 
and is the number of SGs per bin). Although this metric does not properly portray the 
SG spatial density, it serves to illustrate the rapid fall in galaxy counts at small projected 
separations {AR < 150/i~^kpc) due to the fiber colhsion effect. In reality, of course, galaxy 
density {dN/ dr) is a steeply dropping function of separation. Satellite density also falls off 
steeply with increasing velocity separation {AV, Fig. |Dd). 

Figure Ht compares the /raci^e F distribution of HGs and SGs in our group sample to 
that of our entire data set which also includes ungrouped galaxies, i.e. isolated galaxies or 
galaxies near the survey boundary. Relative to the ungrouped objects, HGs are more likely 
to have high /rac-DeF(> 0.9) and slighly less likely to have very low /racZ)eF(< 0.1). SGs in 
groups show the opposite trend: there are fewer high fracDeV {> 0.9) galaxies and more low 
fracDeV (< 0.1) galaxies relative to the ungrouped sample. Neither HGs nor SGs deviate 
noticeably from the parent population for 0.1 < fracDeV < 0.9. Lastly, Fig. Hll compares 
isophotal and moments-based axis ratios of HGs and SGs. We find (in order of increasing 
ellipticity) that HG qmom > SG qmom > HG qiso > SG qiso- We obtain this same order 
using both mean and median qiso and qmom- Using qmom we also find, however, that a larger 
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fraction of SGs exhibit extreme elongation (g < 0.2) than HGs. 

5. Results 

5.1. Analysis Methods and Significance 

To date, most galaxy alignment studies have focused on trend detection rather than 
characterization. We wish to verify and extend these previous analyses by quantifying the 
strength of all alignment signals we find. To that end, we employ the following procedure. 
Within our grouped galaxy catalog, we apply various cuts in different observables (e.g. color, 
AV, AR and fracDeV). The resulting subsamples are binned into histograms to which 
we apply weighted linear least-squares fits. Fit parameters provide measures of both trend 
strength (line slope) and significance (slope uncertainty). We then verify our estimates of the 
significance of the trends we find by comparing to simple analytic estimates. In quantifying 
the magnitude of alignment effects, we are better equipped to isolate the galaxy populations 
responsible for these alignments. 

Our histogram-based analysis affords several advantages. Firstly, linear fits to binned 
data are simple to perform and provide very reliable results with our typical subsample sizes. 
Because of this simple approach, an analytic estimate of the statistical significance of any 
alignment trend is straightforward to calculate. More importantly, by capturing the shape 
of the histogram in a single value (the fractional slope), we are able to concisely investigate 
alignment signal strength in the multidimensional parameter space of galaxy attributes. 
Fortunately, all of the alignment trends we observe are very well described by a linear trend. 
Lastly, the histograms we produce are visually similar to, and thus directly comparable to, 
the results of many previous studies. 

For this analysis, we define a galaxy alignment trend as a systematic variation in the 
linear density of galaxies (the histogram) as a function of either 6, the angular position of 
SGs relative to the HG major axis, or 0, the orientation of the SG relative to the HG-SG 
separation vector ( §3.ip . A fiat (zero slope) histogram implies no alignment signal (isotropy) 
whereas a very steep slope constitutes a strong alignment effect. We use fractional measure- 
ments to better compare alignment trends in galaxy subsamples of different sizes. 

Our trend- fitting procedure operates as follows. First, we define a subsample of grouped 
objects with the desired criteria. We then produce a histogram by binning these objects in 
separation angle. Assuming Poisson errors, each bin i, which contains Ni galaxies, has an 
uncertainty (Tj = v^iV^. Dividing each bin and uncertainty by the bin width yields galaxy 
line density (pi) in each bin. We obtain the best-fit slope {A), its uncertainty {a a), and 
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intercept (B) with a weighted least-squares [x^ minimization) fit to these density values. 
From these values we compute the fractional slope {A/B) and alignment significance {A/aA). 
For convenience, we convert these to percent galaxy density change across the angular interval 
[0,90°], 

AN{%) = 100% (90 °A/B) (9) 

and its uncertainty, 

(TAiv(%) = 100% (90V^/fi) (10) 

By computing the above values in different galaxy subsets, we can quickly identify where in 
parameter space a given effect originates and whether or not variations between subsamples 
are statistically significant. 

To find the source of an alignment effect, we observe how alignment strength and signifi- 
cance vary when the parent population is split into subsamples. For each galaxy attribute we 
investigate (e.g., fracDeV, color), we divide the group sample into different bins according to 
that attribute. If a particular subpopulation is responsible for the observed alignment trend, 
its alignment strength will be diluted by the inclusion of galaxies from other populations. If, 
in dividing the group sample, we improve our selection accuracy (i.e., we create a subsample 
wherein a larger fraction of galaxies produce the alignment trend), we observe an increase 
in alignment strength (steepening of histogram slope). If, however, the resulting subsamples 
all have similar strength to each other or to the parent population, then the alignment effect 
does not depend on that particular attribute used to subdivide them. By repeating this 
procedure we can identify what galaxy subpopulations do and do not exhibit different types 
of alignment. 

In representing individual histograms by a single number (the fractional galaxy density 
change), we are able to expand our histogram-based analysis to multiple parameters simul- 
taneously. We illustrate this technique in Figure O The final result (Fig. divides the 
subsample of highly concentrated host and satellite galaxies {fracDeV > 0.9) into 10 bins 
of radial separation (100/i~^kpc per bin). We separately fit the galaxy density histogram 
and compute the fractional change within each subdivision (Figures [5)d1 - [5)d10) and plot 
these findings as a funtion of radial separation (Fig. [5^). The error bars are in the condensed 
result panel are the slope uncertainties from the individual histogram fits. Applying this 
technique for many different subsamples of observable host and satellite galaxy properties 
(e.g., concentration or color) allows us to identify the more complex parameter dependence 
of certain alignment trends. 

For a given subsample with Nfot galaxies, the following analytic estimate determines the 
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expected alignment significance: 



S 
N 




3 VT^t 



(11) 



where = 45 ° is the half width of the possible range of alignment angles (see Appendix lAl) . 
We find excellent agreement between this estimate and the significance determined explicitly 
by our fitting procedure. 

We now summarize our results for the major axis preference and hostward elongation. 
We compare our results with those from the literature in ^ 



The most noticeable anisotropy we find is a satellite galaxy (SG) distribution which 
favors the host galaxy (HG) projected major axes in projection. For the full group sample 
we observe a — 12.4±0.8% change in galaxy density between major and minor axis at 15a 
significance (black line in Figures IMTTl) . As we will now show, this anisotropy is strongest 
for red, centrally-concentrated SGs of red, centrally-concentrated HGs. Further, the effect 
appears to increase in strength and significance with decreasing physical (3-D) separation. 

Selecting host-satellite pairs by host fracDeV sheds more light on the origin of this effect 
(Fig. ED- While our catalog as a whole exhibits this anisotropy at the 15a level, we find that it 
arises entirely from systems with HG fracDeV > 0.5. More specifically, the signal is strongest 
among systems with HG fracDeV> 0.9 which show a —23.9 ± 1.1% SG density change from 
major to minor axis at 21.2a. Galaxy groups with 0.5 < B.G fracDeV < 0.9 show a weaker 
preference (—6.8 ± 1.6%, 4.3a). At the same time, our fit to this group shows a 3.6-fold 
decrease in significance relative to the total galaxy population. We see no significant trend 
for HG fracDeV< 0.5. 

Splitting the group sample by SG fracDeV, we observe a decreased range of alignment 
strengths (Fig. [7]). The two de Vaucouleurs profile-dominated SG subsamples (SG fracDeV> 
0.9 and 0.9 > SG fracDeV > 0.5) show enhanced alignment strength (—17.2 ± 1.3% and 
— 15.9 ± 1.8%) relative to the total population and are consistent with each other. The two 
exponential profile-dominated SG subsamples (0.1 < SG fracDeV < 0.5 and SG fracDeV < 
0.1) are also consistent with each other (changes of — 7.1±1.8% and — 6.2±1.7%, respectively) 
but shallower than the combined group sample. 

Breaking the sample down by HG and SG color shows a similar picture. We observe 
a very strong effect between red HG - red SG pairs (—26.2 ± 1.3%, 19.5cr, Fig. [8]). Blue 



5.2. Major Axis Preference 
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SGs around red HGs show a less significant alignment trend (—14.3 ± 1.5%, 9.6cr). We find 
no evidence of anisotropy among either red or blue SGs of blue HGs. These latter two 
subsamples are consistent with each other (—1.5 ± 1.8% and —2.6 ± 2.2% decreases) and 
with isotropy. We deduce that red and centrally-concentrated (i.e. elliptical) HGs and SGs 
are almost exclusively responsible for the major axis alignment preference. 

We observe this major axis preference across a wide range of HG-SG velocity separation 
(Fig. E>- Although it is most significant for \AV\ < 300km/s (-14.8 ± 1.0%, 14.3cr), we 
observe an equally strong (14.8% decrease) trend for 300km/s < |AV^| < 600km/s albeit 
at somewhat reduced significance {8.3a) due to the smaller subsample size. The subsamples 
at greater velocity separations {AV > 600km/s) are consistent with each other and with 
isotropy. We find a similar effect in projected radial separation (Fig. [TU]) . Alignment is 
strongest (-20.4 ± 2.0%, 10a) for the closest SGs {AR < 250/i-^kpc) but is still detected 
at nearly 6a to the group size limit of 1 Mpc. Taken together, these velocity and radial 
separation trends indicate an increasing spatial anisotropy in favor of projected HG major 
axis with decreasing 3-D separation. 

We also investigate anisotropy as a function of (spectroscopic) galaxy group size (Fig. 
[TT]) . We find that alignment strength increases steadily with increasing group size. Among 
SGs from the smallest groups (1 < Nsg < 5; 95186 objects) we observe a —9.6 ±1.1% galaxy 
density change towards the minor axis at 9a significance. In the largest groups {Nsg > 20; 
16122 objects) we find a —23.5 ± 2.4% change (9. 60"). Interestingly, the latter subsample 
exhibits the strongest alignment and highest significance despite having the fewest galaxies. 
From this we deduce that major axis preference exists for all group sizes but is strongest for 
the largest galaxy groups. 

In Figure [T2] we explore the interdependence of physical separation and light profile 
type on the observed alignment trend. We limit velocity separation to |A1^| < 500km/s and 
explore the dependence of alignment strength on projected radial separation {AR). We find 
that when both HG and SG have high /racDe F (especially fracDeV > 0.9, top-right panel), 
the anisotropy strength (fractional slope) increases with decreasing AR. We conclude that 
the dependence of anisotropy on physical separation occurs largely among the same red, 
centrally-concentrated galaxy pairs described above. We further discuss the significance of 
these results in §7.11 
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5.3. Hostward Elongation 

We find that satellites tend to elongate towards their host galaxies. For the full group 
sample, we observe a —5.7 ± 0.8% change (at 7a) between the frequency of satellites that 
exhibit perfect radial (hostward) alignment relative to those that exhibit tangential alignment 
(black histogram in Figures [TSlfTS]) . Of the galaxy properties we examined, the degree of 
alignment is most strongly dependent on SG fracDeV (concentration. Fig. [T3l) . SGs with 
fracDeV > 0.9 (52,817 objects, 29% of sample) exhibit a —13.9 ± 1.4% density change 
favoring direct hostward alignment at 9.9cr. None of the other SG /rac-De F subsamples 
exhibit significant alignment. 

We present hostward elongation (</>) for different HG profile types in Figure [TH Satellites 
of high-fracDeV EGs (HG fracDeV > 0.9, 94,716 objects, 52% of sample) show a -5.8 ± 
1.1% change in density at 5.3cr. Satellites of HGs with 0.5 < fmcDeV < 0.9 show only 
a marginal (3.9cr) —6.2 ± 1.6% change. The slopes of the trends seen in all four bins are 
statistically consistent with the total group sample, suggesting that hostward elongation is 
largely independent of HG fracDe V. 

Galaxy colors (Fig. [T5|) have some effect on alignment strength. Red SGs around red 
HGs (56,421 objects, 31% of sample) show a —9.8 ± 1.4% density change at 7a. Red SGs 
around blue HGs exhibit a slightly steeper trend of —11.3 ± 1.9% at somewhat lower signif- 
icance (5.8(j). These two subsamples are statistically consistent, providing further evidence 
that hostward elongation is principally dependent on SG properties. We do not observe a 
statistically significant alignment trend among blue SGs of either blue or red HGs. 

There is no evidence that the degree of hostward elongation depends on velocity sep- 
aration (see Fig. [T6|l . All four subsamples are essentially statistically consistent with the 
full group sample. In projected radial separation (Fig. [17]), we find marginal {5a) evidence 
for hostward elongation in both the largest AR and smallest AR bins. At large separation 
(750/i~^kpc < AR < 1 Mpc) we find a —7.2 ±1.4% density change at 5.3cr. In the inner- 
most radius bin [AR < 250/i~^kpc), we find a — 9.7± 1.9% density change at 5. la. These 
results suggest that 3-D physical separation plays a mininal role in the degree of hostward 
elongation. 

Lastly, we examine the influence of group size on hostward elongation (Fig. [T5]) . All four 
group size subsamples are consistent with the whole group sample, suggesting that hostward 
elongation is independent of group size. 
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6. Caveats and Complications 

6.1. Discrepant Position Angles 

Our greatest cause for concern is the accuracy of measured galaxy position angles (PAs). 
During our analysis, we identified systematic PA discrepancies between model- and isophote- 
based position angles (PAdcv and PAjso) that critically affect our results. Some of these dif- 
ferences, primarily those due to isophotal twisting, are physically genuine. However, in many 
cases, isophotal contamination by artifacts and nearby objects appears to be the culprit. We 
suspect that these variations may be largely responsible for the general disagreement be- 
tween many of the results on hostward elongation presented in the literature. We detail our 
findings in the remainder of this section. 

To investigate this PA discrepancy, we compare alignment signals of identical galaxies 
using PAiso and PAjey. In order to illustrate the discrepancy, we do not enforce agreement 
between the two measured PA values. We do, however, continue to require that AR > 4rp. 
We first examine location angle [6) which requires accurate host galaxy (HG) PA. Comparing 
our PAiso- and PAjev-based results (Fig. [2D]), we observe few significant differences. 

The situation changes dramatically for the hostward elongation case (Fig. [21]). Isophotal 
and model PAs produce markedly different results in high fracDeV SGs (rightmost column), 
particularly those of high fracDeV B.Gs {fracDeV > 0.9, top-right panel). The strong (and 
high-significance) hostward alignment we see with PAiso at small projected radii {AR < 
100 h^^ kpc) is absent from PAdcv-based measurements of the same satellite galaxies. Figure 
[22] compares side-by-side the individual histograms of this highly discrepant subsample before 
and after applying the APA < 15 ° criterion. After rejecting SGs with APA > 15 ° (30% of 
objects), both the alignment trend and discrepancy are significantly reduced. These results 
indicate that the apparent degree of hostward elongation is quite sensitive to PA selection. 
How best to deal with this depends on the cause of the discrepancy, which we now investigate. 

Neither PA measure is perfect. The galaxy models (exponential and de Vaucouleurs) are 
fitted to each galaxy across a range of radius. These fits account for seeing effects (^ and 
may interpolate across artifacts and substructures as necessary, which reduces the potential 
for errors due to contamination. There are, however, several known yet unsolved problems 
within the model-fitting routine. Most importantly, the galaxy model PA distribution is 
significantly enhanced between 60 ° and 145 ° (nearly East- West in orientation) rather than 
uniformly distributed in angle (Fig. [23]). In addition, th e derived galaxy scale radii (an d 



thus the derived axis ratios) exhibit some discretization (lAdelman-McCarthy et al.[[2008[ ). 
Isophotal PAs, by contrast, are determined by a best-fit ellipse to 25 mag / sq. arcsec surface 
brightness around each galaxy. Unlike the matched galaxy models, the isophotal ellipse fit 
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is sensitive to a narrow and faint surface brightness region of each galaxy. As such, the 
isophotal eUipses (and thus their measured PAs) are more susceptible to contamination by 
nearby objects and image artifacts. 

To explore the cause of these discrepant PA measurements in individual systems, we 
obtained SDSS images of all high-fracDeV {fmcDeV > 0.9) SGs hosted by high-fracDeV 
HGs with discrepant hostward elongation signals. We visually compared the SDSS-reported 
isophotal ellipses and de Vaucouleurs PAs and scale radii with the galaxy images. We 
found that the great majority of discrepant PAs could be attributed to two causes: intrin- 
sic variations in the PAs of isophotes with radius (~ 70%) and PA misestimation due to 
contamination (~ 30%). Contamination was caused (in roughly equal numbers) either by a 
nearby bright object or by a bright knot somewhere within the galaxy. 

Figure [21] shows an example of an isolated galaxy for which the model and isophotal 
PA differ by nearly 90 ° due to intrinsic isophotal twisting. Generally, the de Vaucouleurs 
and exponential galaxy profiles fit to the brighter galaxy interior while the 25 magnitudes 
/ arcsec^ isophote traces the outer edge of each galaxy. This is particularly noticeable 
in this example, where the de Vaucouleurs model fits primarily the central bar or bulge 
region, while the isophotal ellipse traces the outer edge of the spiral arms. In this fashion, 
isophotal twisting can lead to highly discrepant PA measurements. Given that the physical 
mechanisms which give rise to galaxy alignments are not well understood, it is not clear 
which PA determination is most relevant. This motivated our decision to exclude all galaxies 
with highly discrepant (APA > 15 °) isophotal and de Vaucouleurs PAs from this analysis. 
Further study is required before these interesting objects can be properly incorporated into 
an alignment analysis. 

Figure [22] shows illustrative examples of isophotal contamination. In a number of cases, 
failure to separate the light from other objects results in an isophotal PA which does not 
physically reflect the system in question. Although not as prevalent as isophotal twisting, 
such cases are responsible for a significant fraction of discrepant PA measurements. 

Importantly, we note that we do not expect random (i.e., isotropic) PA contamination 
for galaxies in groups and clusters. In such environments, galaxy spatial density is a rapidly 
declining function of separation from the center (AR). As a result, a given group or cluster 
member is most likely to have an optical companion radially interior to its location. We thus 
expect that isophotal contamination due to chance juxtaposition of galaxies will preferentially 
skew isophotal PA in the radial (hostward) direction. The degree to which this occurs will 
depend strongly on the how rapidly galaxy number density falls off with increasing separation 
from the cluster center. 
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From this investigation we conclude that care is required to ensure that only accurate 
PA measurements are included in alignment analysis. This is particularly important for con- 
tamination cases where the best-fit isophotal ellipse is highly elongated, giving the illusion of 
high accuracy. Elimination of discrepant cases has a profound effect on the measured host- 
ward alignment signal. We believe the use of different PA measurements may be responsible 
for several conflicting previous results (see discussion in §7.21) . 



6.2. Spectroscopic Survey Limitations 



In some cases, the hmitations of the SDSS spectroscopic sample restrict our ability to 
properly characterize alignment. Fiber collisions make it difficult to observe galaxies with 
small projected separations. On SDSS plates, the spectroscopic fibers (3" each) can be set 
no closer than 55" apart. Although this limit is mitigated somewhat by plate tiling, many 
of the closest neighbors in projection will be absent. In addition, the finite (640) number 
of fibers allowed per plate ultimately limits the fraction of galaxies that can be observed 
in dense fields. As a result, the completeness of a galaxy group generally decreases with 
increasing group size. The inability to fully sample the innermost regions of groups and 
clusters could systematically affect the measured alignment signal. Of importance to our 
exploration of the Holmberg Effect in the SDSS spectroscopic sample, some bright nearby 
neighbors may be absent because of these effects. If such omissions were common, we would 
incorrectly identify many objects as isolated galaxies. Whether or not the robust elimination 
of interlopers afforded by spectroscopic redshifts outweighs the potential statistical benefits 
of the larger photometric sample for detecting alignment trends is a difficult question which 
we do not attempt to answer here. It is likely, however, that exclusive use of galaxies from 
the spectroscopic sample will not provide the most accurate picture of galaxy alignment at 
separations approaching the fiber collision limit. 

The magnitude limit of the spectroscopic survey may also affect our results. It was 
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the apparent shape of a cluster (isopleths) may vary as a function of the galaxy flux limit, 
so that the distribution of the brightest members may not be representative of the overall 
population. A similar effect due to either the spectroscopic magnitude limit or fiber collisions 
is possible with these data. 

The survey magnit ude limit al s o like ly hinders our ability to properly reproduce the orig- 
inal sel ection criteria of iHolmbera fjl969l ) and robustly assess the existence of the Holmberg 
Effect. iHolmberg] (119691 ) mandated a host-satellite mass ratio of at least 25 (corresponding 
to a flux ratio of roughly 3.5 magnitudes). With the limiting SDSS survey magnitude of 
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17.77, only host galaxies of r < 14 are bright enough to have spectroscopically-identified 
"Holmberg-like" companions. Of the 572,495 galaxies in our dataset, only 2,775 are this 
bright. We find only 1178 such HGs in our group catalog. 

To make matters w orse, scattered light in the SPS S telescope increases sky brightness 
around bright galaxies ( Adelman-McCarthy et ahlboosl ). This effect is believed to conceal 
fainter nearby galaxies, resulting in an observed dearth of faint galaxies near bright galaxies 
( lAdelman-McCarthy et al.ll2008l ). Several recent studies of weak lensing and SG anisotropy 
mask objects at very close separations to avoid these potential sele ction biases (for example, 
obje cts in the innermo st 30 kpc and 35 kpc are excluded in iMandelbaum et al.ll2005 
and lBailin et al.ll2007l . respectively). In particular, this seriously complicates application of 
isolation criteria. 



6.3. Isolation 



Well-selected isolation criteria are of equal importance to major axis preference and to 
hostward elongation, both of which have now been obs erved to varying degrees in multiple in- 
dependent investigation s. Several of these studies (e.g., Brainerdll2005 : Agustsson fc Brainerd 
20061 : lYang et al.ll2006l ) employ aggressive selection criteria in order to extract a sample of 
isolated objects. Among these isolated galaxy systems, they observe strong SG distribu- 
tion anisotropy favoring the HG major axis. The aforementioned studies additionally report 
near-perfect SG isotropy around blue hosts (i.e., no sign of the Holmberg Effect). 



Recently, iBailin et al.l (120081 ) performed a systematic evaluation of the various selection 
criteria used in recent alignment studies. They conclude, based on examination of mock 
catalogs and simulations, that most previous efforts to select only isolated galaxy systems 
were ineffective. Objects people believed were truly isolated systems were statistically shown 
to include many group and cluster members. In reality, the SDSS spectroscopic survey alone 
does not, in general, include all potential neighbor galaxies, which makes application of 
isolation criteria very difficult. Based on our satellite anisotropy results ( §5.21) and the HG 
and SG properties most closely associated with the alignment signal, we conclude that the 
alignment signal is largely due to the contribution of group and cluster members despite the 
ostensibly isolated nature of our group samples. 
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7. Previous Detections and Comparisons 
7.1. Satellite Galaxy Anisotropy 

The tendency for satellite galaxies to have anisotropic distributions favoring t he major 
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Bailin et al. 2008 ) and with decreasing radial separation fjYang et al. 20061 ) have been inde- 



pendently observed. Our results are consistent with these other recent findings in both in 
strength and nature with these other recent findings. 

We build on these results by demonstrating the monotonic dependence of major axis 
preference on HG and SG light profile concentration, physical separation (both AV and 
AR) and group size. In §5.21 we show that red, highly concentrated SGs of red, highly 
concentrated HGs preferentially reside near the projected major axis of their hosts. We also 
show that this effect increases in strength with decreasing physical separation. Lastly, but 
equally importantly, we observe that the alignment strength increases with increasing group 
size. 

Galaxy clusters are known to produce a similar effect. Many past investigations have 
discovered elongated galaxy cluster is opleths which align strongly with the rn ajor axis of the 
brightest cluster galaxy (BCG) (e.g. ICarter fc Metcalfelll98ol : lBinggelilll982h . Indeed, BCG 
- cluster alignment is sufficiently strong that the BCG PA can be used as a proxy for the PA 
of the cluster galaxy distribution. 

As would be expected in cluster environments, the galaxies which give rise to our mea- 
sured anisotropy are mainly red, highly concentrated ellipticals. In the literature we find 
good agreement between our alignment fi ndings and those from de ep c ampaigns of in divid- 
ual galaxy clusters (e.g., the histograms of ICarter &: Metcalfelll980l and lBinggelilll982l ) . This 
further supports our suspicion that group and cluster members produce the majority, if not 
all, of the observed anisotropy signal. 



7.2. Hostward Elongation 

We observe hostward elongation in our SG sample ( §5.31) at the few percent level (mea- 
sured in fractional galaxy density change), in fair agreement with numerous previous studies. 
We find that the alignment signal is strongest among high fracDeV {fracDeV > 0.9) SGs 
irrespective of host type. We also observe an enhancement among red SGs (independent 
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of HG color). Most importantly, we report a connection between the type of PA used and 
the apparent alignment strength which may explain the variation among several previous 
results. 

For some time, authors have suspected that intrinsic galaxy a li gnment effe c ts ma y act 
as a contaminant in weak lensing surveys (jCroft fc Metzlerl l2000l ). iLee fc PenI (120011 ) pre- 
dicted that intrinsic alignmen t s wou ld dominate the weak lensing signal for spiral galaxies in 
SDSS. iBernstein fc Norberg I (120021 ) found no significant tangential elon gation in a study o f 
1819 galaxies from the 2dFGRS, limiting the contamination to ~ 20%. iHirata et al.l (120041 ) 
found no evidence fo r intrinsic alignment of s at ellite galaxies measured w ith the adaptive 
moments technique ( Bernstein fc JarvisI I2OO2] ). iMandelbaum et aL ( 2005 ) also concluded 
that contamination due to intrinsic galaxy alignment is minimal. 

Using our entire SO sample, we find that the frequency of SGs oriented radially towards 
their hosts is 6% greater than that of SGs oriented perpendicular to their hosts ( §5.31) . This 
anisotropy is strongest among high-/racDeF SGs, for which the fractional change steepens 
to —1 4%. These findin g s are generally consistent with those mentioned above. On the other 
hand, jPereira fc KuhnI (120051 ) report a much stronger radial alignment in SDSS data with 
a sample of 85 X-ray s electe d clusters. They observe a —25.5% change at Aa significance. 
Agustsson &: BrainerdI (120061 ) also find a much stronger radial elongation effect. They find 
that SGs located at small projected separation [AR < 100 kpc) from relatively iso- 
lated hosts exhibit a ~ —20% density change between radial and tangential orientation. 



Faltenbacher et al.l (120071 ) reports a similarly strong effect for red satellites at low projected 
separation. 

The last three studies, mentioned above, relied exclusively on isophotal PA measure- 
ments, which, we believe, are suspect. We demonstrated the impact of PA type (isophotal 
or model-based) on the measured hostward elongation signal in §6.11 By varying our PA 
selection, we can recover results from a wide range of recent investigations. We conclude 
specifically that relying solely on SDSS isophotal PA produces a stronger radial (hostward) 
elongation signal than does the exclusive use of matched models or adaptive moments. We 
present our interpretations and further discussion in ^ 



8. In Search of the Holmberg Effect 



HolmbergI (119691 ) observed that the distribution of nearby (projected radial separation 
< 40 kpc) satellite galaxies (SGs) around inclined, isolated spiral host galaxies (HGs) 
was heavily concentrated near the HG projected minor axis. In fact, after statistically 
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subtracting the expected isotropiccally-distributed field galaxies, Holmberg found effectively 
no SGs within 30° of the projected HG major axis {6 < 30°). Since the host galaxies were 
inclined spirals, this suggested that SGs of such hosts have primarily polar orbits. 



Holmberd (119691 ) selected isolated systems for several reasons. First, such systems are 
dynamically simple in comparison to larger groups. In addition, stringent isolation crite- 
ria could be used to decrease t he frequen c y of i nterlopers in the SG sample. To evaluate 



a galaxy's degree of isolation, iHolmberg] (119691 ) first defined a circular survey region of 



40/i~^kpc projected radius around each HG. Next, he estimated the (stellar) mass of each 
HG and SG companion from luminosity and morphological type. Isolation required that all 
SGs within twice the survey radius (80/i~^kpc) have less than 1/5 the mass of the HG and 
all SGs within the survey radius {40h~^kpc) have less than 1/25 the HG mass. Systems 
which passed those criteria and were sufficiently inclined {q < 0.53) were included in the 
analysis. 

In practice, exactly replicating these conditions using the SDSS spectroscopic galaxy 
sample is essentially impossible, primarily because of spatial sampling constraints (i.e., "fiber 
collisions," see §6.2p . Insufficient spatial sampling causes multiple problems. Firstly, the very 
small typical separations of Holmberg's original SG sample (40 kpc is roughly 4x the HG 
semi-major axis length) greatly intensifies the incompleteness of the sample due to the fiber 
collision effect. After removing likely blends and misidentifications (SGs with AR < Vp or 
AR < isoji) we find only 132 SGs around 130 inclined {qiso < 0.53 and qmom < 0.53) HGs in 
our original group sample (i.e. without any isolation criteria applied) with AR < 40 kpc. 
Without the ability to identify all nearby galaxies of comparable size and brightness, it is 
impossible to properly apply Holmberg's isolation criteria to our galaxy sample. 



8.1. An Isolated Sample 

We first attempt to apply Holmberg's selection criterion (approximately) to our group 
sample in order to demonstrate the difficulties with reproducing his analysis precisely with 
the SDSS spectroscopic sample. Unable to exactly reproduce his galaxy mass estimates, 
we apply the same cuts using bolometric luminosity (calculated with the LRT codes). The 
modified group selection procedure incorporates the following criteria: 

{Ltoiho > 25 {L,,i)sG for AR < 40 kpc (12) 
(Lboihc > 5 iLboi)sG for AR < 80 h-' kpc. (13) 

All hosts which fail these tests are removed from the HG list. The remaining groups comprise 
our 'isolated' sample. We identify inclined HGs by axis ratio, requiring that qmom < 0.53 and 
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Qiso < 0.53. Unfortu nately, this proce dure yields only 2 SGs of 2 HGs within the neighbor 



volume employed by iHolmbera (jl969l ). 



8.2. More Distant Companions: a Marginal Effect 



Unable to replicate the analysis of iHolmbergj (Il969l ) with SDSS spectroscopic data, we 



instead proceed with an examination of more distant SG companions. We select SGs of 
inclined isolated systems as above but relax the radial separation constraint, including SGs 
as distant as 500 kpc. In the spirit of selecting relatively isolated HGs, we exclude SGs 
from systems with more than 25 detected satellites within 1 Mpc. From among these, we 
select the SGs with 15 x SG L^oi < B.G L^oi and SG fracDeV < 0.10. These last few criteria 
specifically pick out smaller SGs while avoiding clusters and big groups. 

Applying the above criteria, we observe a preference for the minor axis among the 
remaining SGs (Fig. UM . Among SGs with AR < 500 kpc and AV^ < 500km/s (51 
galaxies), we find few galaxies (4 galaxies per 15° bin) near the HG major axis [6 < 30°). 
The density increases to 9 galaxies per bin for 30 ° < ^ < 60 ° and to 12 and 13 galaxies 
per bin for 60° < 6* < 75° and 75° < 6' < 90°, respectively. This roughly 3-fold increase 
indicates anisotropy at 3cr significance. For SGs with AR < 300 kpc and AV < 300 km/ s 
(27 galaxies), we find a low, constant galaxy density (2 galaxies per 15° bin) for 6* < 45°. 
For 6' > 45°, we find that galaxy density rises steeply, such that galaxies with 75° < < 
90° are 4.5 times as numerous (9 galaxies) as those with 6 < 45°. For the latter case, 
although the fractional slope is quite large, the dwindling galaxy counts ultimately reduce 
the significance of the deviation from isotropy to approximately 2.5a. Keep in mind that 
rather than including all of the satellites in the galaxy systems, we have selectively included 
and excluded SGs with different properties to obtain this behavior. We caution that that 
our results depend quite sensitively on the specific selection criteria we use. With so few 
objects, further subdivision is not practical. 

In summary, we observe that a subset of SGs with certain properties favor the HG minor 
axis within 500 h^^ kpc projected radius at marginal significance. The strength of this trend 
may increase with decreasing AR, but the accompanying decrease in the total number of 
SGs decreases the significance of this result. Whether or not this apparent anisotropy and 
preference for the HG minor axis is due entirely to selection effects remains to be seen. 
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8.3. Difficulties and Discussion 



Although we observe a tren d that resembles the Holmberg Effect, we are un a ble to 
confirm or refute the findings of iHolmbergl (119691 ) in the present work. iHolmberg fjl969l ) 
enforced isolation criteria using a mass cut. In practice, his mass determination is difficult 
to reproduce. As a substitute, we impose the same cuts in bolometric luminosity {Lboi 
computed with LRT). We also use SDSS fracDeV and K-corrected galaxy colors in lieu of 
confirmed morphological type. 

A more irn portant difference is the spatial extent of satellite survey regions. The 
Holmberg] (119691 ) sample is restricted to faint nearby (typically within 4 galaxy radii) com- 
panions. Very few such galaxies even exi st in the SDSS spectro scopic sarn p le (^on ly find 2 
after imposing criteria of Holmberg! Il969l ). At close separations, Holmberg (1969) requires 
a host-satellite mass ratio of at least 25. In practice, many such galaxies will be near the 
magnitude limit of the spectroscopic survey. Those that are sufficiently bright may be un- 
dersampled because of fiber collisions (see 



Zaritsky et al.1 (Il997br) detect a Homberg-like SG anisotropy using their own indepen- 



dent sample (jZaritsky et al.lll997al ) at larger radii (> 200 kpc). They find inconsistency with 
isotropy at the 99% confidence level, reporting that 72 of their 115 SGs (63%) reside at 
^ > 45°. This corresponds to a 217% fractional increase in the number of SGs aligned with 
the minor axis relative to the major axis, assuming the underlying distribution is linear in 

e. 



We find an alignment signal resembling that of iHolmberg] (ll969f ). albeit at larger pro- 
jected s eparations. In magnitu de and radial extent, our results are generally consistent with 
those of I Zaritsky et al.l (ll997bl ). although at marginal significance. We further caution that 
our findings are highly sensitive to selection criteria and may be the result of selection effects. 



9. Discussion and Conclusions 

We conclusively observe an anisotropy in the distribution of satellite galaxies (SGs) in 
systems with red, highly concentrated (chiefly elliptical) host galaxies (HGs). Within these 
systems, the alignment trend strength is further enhanced by nearly a factor of two for red, 
high-concentration SGs (also mainly ellipticals). We find conclusive evidence of increasing 
alignment strength with decreasing physical separation (both in projected radius AR and 
velocity separation AV^), noting that the fractional anisotropy among the closest satellites 
[AR < 250 h"^ kpc) is 40% greater than that seen among more distant (500 kpc < AR < 
1000 kpc) companions. In addition, the strength of the alignment trend increases with 
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spectroscopically-identified group size (see §5.21) . 



Although these results agree well with several other recent studies (see §7. II) . we offer a 
slightly different interpretation. Based on the nature of the galaxies which give rise to this 
effect, we suspect that members of large groups and clusters contribute significantly to, or 
even dominate, this alignment trend. The well-known alignment of the members of galaxy 
clusters with the orientation of the brightest member has a c omparable rnagnitu de and sign. 
A recent thorough inspection of popular selection criteria (IBailin et al.l |2008| ) determined 
that extremely stringent criteria are required to generate a sample of truly isolated hosts 
and satellites. They find that with less stringent criteria, group and cluster members will be 
ubiquitous in ostensibly isolated SDSS galaxy samples. 

The cluster origin offers a natural explanation for this alignment effect. Cluster shape 
is defined by the spatial distribution of member galaxies. These members form in the high- 
density filaments, creating an anisotropy in the distribution of satellites during the formation 
process. Although dynamical processes may alter the distribution somewhat as time passes, 
numerical work indicates that the cluster shape w ill r etain the memory of th e initial member 
galaxy distribution (see e.g., iKnebe et al.l |2004| and iLibeskind et al.ll2005l ). The brightest 
cluster galaxies (BCGs) could acquire their alignment with this anisotropic distribution 
of member satellites through tidal forces exerted by the large-scale structure during their 
formation. Alternatively, mergers and accretion e vents drawin g frorn the already-anisotropic 
SG distribution could create this alignment. See iBailin et al.l (120081 ) for further details and 
discussion of this scenario. 

We also find significant evidence of hostward elongation among satellite galaxies (SGs). 
Although this effect has been observed before, we are able to reconcile many previously 
conflicting results. Among the attributes we investigate, hostward elongation depends most 
strongly on SG concentration and color. We find effectively no significant direct dependence 
on HG properties. Ultimately we conclude that the hostward elongation is real, but not 
nearly as strong an effect as the major axis anisotropy mentioned above. 

More importantly, we find that the detection of hostward elongation is a strong function 
of the method used to measure and define position angle of the galaxy. We find that, at 
low projected separation (Ai? < 100/i~^kpc), the radial alignment signal calculated with 
isophotal PA sharply increases in strength. We observe no such increase among the same 
galaxies using de Vaucouleurs PA. We further find that excluding objects whose PA measure- 
ments differ_by_jnoreth^;^ This explains why some a uthors 



Pereira fc KuhnI |2005| . lAgustsson fc BrainerdI |2006| . & iFaltenbacher et al.l 120071 ) , who 



(e.g. 

relied solely on isophotal PA for galaxy orientati on, observed a stronger rad i al alignment 



effect than those who used galaxy moments (e.g., iBernstein fc Norberg 1120021 . iHirata et al. 
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20041 . & iMandelbaum et al.l |2005| ) . Given that PA determination plays such a crucial role 
in galaxy alignment studies, we stress the prevalence of these discrepant PAs and their 
significant potential impact. 

At the same time, this observed discrepancy raises new questions. Host galaxies (HGs) 
in our group sample appear immune to this phenomenon ( §6.ip . The rate at which SGs 
are affected by this appears to increase dramatically at small projected separations from 
the HG. Observing the anomalous objects by eye, we find that most exhibit strong isophotal 
twisting (e.g. barred spiral galaxies), although there is a significant (nearly 1/3) contribution 
from isophotal contamination. As outlined previously, isophotal contamination may induce a 
radial alignment bias if spatial galaxy density is a decreasing function of projected separation. 

For cases of severe isophotal twisting, position angle is an ill-defined quantity. Since the 
mechanisms of galaxy production and alignment within their haloes is still highly uncertain, 
it is not clear how to best incorporate galaxies with extreme isophotal twisting into an 
alignment investigation such as this. For the same reason, we ar e unsure how to interpret 



a radial elongation tendency. Some numerical simulations (e.g. iFaltenbacher et al.l l2008l ) 
show a strong tidally-induced radial alignment effect in dark matter haloes. Whether or not 
such orientations would also apply to the luminous regions of galaxies within these haloes is 
unclear. 

Regardless of the cause, measuring the intrinsic radial elongation signal of galaxy sys- 
tems is an important step in the applicati on of weak lensing to precision cosmology. In 



gener al agreement with other authors (e.g. iBernstein fc Norberg 1120021 : IMandelbaum et al. 



20051 ). we believe that a radial elongation signal with the magnitude we have found in SDSS 
spectroscopic galaxy sample will minimally influence galaxy-galaxy lensing at the levels of 
precision currently required. Further investigation, however, will be required to determine 
whether or not this hostward elongation is observed among fainter galaxy populations (i.e., 
in the SDSS photometric galaxy sample). 

A third alignment trend, the Holmberg Effect, remains improperly tested. The mag- 
nitude and spatial sampling limitations of the SDSS spectroscopic survey combined with 
systematic effects in the vicinity of bright galaxies effectively preclude implementation of 
Holmberg's original selection criteria. Instead, we explore the distribution of SGs at larger 
projected separation. We flnd that, although the whole SG population around inclined, 
isolated spirals shows no detectable anisotropy, speciflc SG subsets do. SGs from smaller 
systems {Nsg < 25) with fracDeV < 0.15 and L^oi < 1/15 that of their host overwhelmingly 
prefer the HG minor axis, albeit at marginal (3.4(t) signiflcance owing to the relatively small 
number of objects that satisfy these criteria. 
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Through improved convergence of independent observational results, we have uncovered 
and characterized several different galaxy alignment trends. However, further refinement is 
required before high-accuracy alignment prescriptions are available for weak lensing and 
simulations of galaxy formation and evolution. Importantly, much of this remaining work 
pertains to identification and removal of systematic effects such as position angle uncertainty 
that currently prevent a more thorough understanding of intrinsic galaxy alignment effects. 
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A. Expected Signal-to-Noise 

Here we derive the expected signal-to- noise ratio (S/N) of a trend detection using our histogram- 
fitting method ( ^5.ip . Here S/N is the square root of the Ax^ improvement of a weighted hnear 
least-squares fit relative to a constant value (mean). Larger Ax^ values constitute stronger evidence 
for anisotropy. We define: 

S/N = v^A^. (Al) 

In our linear model, the histogram has n^, bins in the angle 6. Each bin i contains 

Ni = A9ide + N (A2) 

galaxies, where N = Ntot/nh is the mean number of galaxies per bin, dO is the bin width, Ntot = 
Yl'i=i total number of galaxies in the sample, and A is the best-fit slope of a linear fit to 

the histogram bins. The constant model assumes that Ni = N. 

The improvement of a linear least-squares fit over the constant value is the squared difference 
between the two models divided by the individual bin counts, summed over all histogram bins. The 
N terms in the linear and constant models cancel, leaving: 

Ax' = ±^-^^. (A3) 

1=1 * 

To simplify the calculation, we assume that the Ni are roughly equal and that Ni ~ N. 
Removing bin-independent terms from the summation, we have: 

^x' = =def20fde. (A4) 

i=l 

To simplify the calculation, we define AO to be half the histogram width (AO = 45 ° in this 
case). We can relate this to bin width (d9) by noting that 2A6 = n^dO. Substituting this outside 
the summation, we have: 

Ax2 = ^f^y02^0. (A5) 

N 

We convert this to an integral in the limit n;, ^ co, leaving 

Ae 

, 2 A^ 2A0 f . A^ 2Ae 2 , „ A^ AA9^ AA^AO^ , , , 

Ay^ = = / e^dO = ■= A9^ = = = . (A6 

N m J N rib 3 N 3nb 3Ntot 

Finally we have, 

S ^— ^ [i AAO^ 
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Fig. 1. — Characteristics of SDSS galaxy parameter fracDeV. Panel (a) shows the fracDeV 
distribution of red and blue galaxies separately as fractions of the total sample (572495 
objects total). The red and blue curves represent red and blue galaxies, respectively. The 
green curve is the total of all galaxies. All galaxies in the initial sample are represented. 
Panel (b) presents the full range of values for inverse concentration (Petrosian R^o/Rgo). 
Different values of fracDeV are isolated to illustrate the distinct groups that form the overall 
curve. The sharp cutoff towards low inverse concentration is partially caused by the effects 
of seeing. Panel (c) shows the strong correlation between fracDeV and inverse concentration 
among galaxies of < fracDeV < 1. 
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Fig. 2. — Galaxy sample col or propertie s. We adopt the locus {g — r) = 0.78 — 
0.0325 (Mr — blogh + 19) of iBailin et al.l (120071 ) (shown in red) as the color division be- 
tween blue cloud and red sequence galaxies. Panel (a): all galaxies in our sample (panels b, 
c, & d combined). Panel (b): pure de Vaucouleurs galaxies {fracDeV = 1). Panel (c): pure 
exponential galaxies {fracDeV = 0). Panel (d): combined-profile objects (0 < fracDeV < 1). 
The gray scale is normalized seperately in each panel. 
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Fig. 3. — Left: Neighbor volume cylinder diagram. Satellite galaxies have AR < 1 Mpc 
and AV^ < 1200 km/s relative to their host galaxies. Right: schematic diagram illustrating 
the angles 6 and we use to characterize the relative positions and orientations of host and 
satellite galaxies. The angle 9, also called location angle, measures the angular position of 
a satellite galaxy relative to the direction of the host galaxy major axis. measures the 
orientation of the satellite major axis relative to "hostward" (radial). 
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Fig. 4. — Properties of the galaxy group sample. Panel (a): satellite galaxy (SG) counts as a 
function of projected radius [AR). Panel (b): SG density as a function of velocity separation 
(AV). Panel (c): host galaxy (HG) and SG /racL'e F relative to the ungrouped data set. 
Panel (d): isophotal and moments axis ratios {qiso, Qmom) for both HG and SG populations. 
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Fig. 5. — Illustration of our histogram condensation analysis technique. The original sub- 
sample, consisting of hosts and satellites with high concentration [fracDeV > 0.9), is further 
subdivided into 10 bins of radial host-satellite separation (100/i~^kpc per bin, panels bl- 
blO). We fit a linear model to this histogram within each subdivision and compute the 
fractional change and its uncertainty. These fractional changes are then reproduced in the 
panel (a) as a function of radial separation. The error bars are the slope uncertainties from 
the fitting procedure. 
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Host Major Axis Preference: host galaxy concentration 







HG frocDeV 




AN [0:90] 






|A|/<7, 






Total : 


162954 


-12.4% 


± 


0.8% 


15.36 




HG 


fracDeV > 0.9 : 


73395 


-23.9% 


± 


1.1% 


21.16 


0.5 


< HG 


frocDeV < 0.9 : 


44319 


-6.8% 


± 


1.6% 


4.26 


0.1 


< HG 


fracDeV < 0.5 : 


27722 


3.0% 


± 


2.1% 


1.44 




HG 


fracDeV < 0.1 : 


17518 


4.8% 


± 


2.7% 


1.80 



CD 




20 40 60 80 

Location Angle (6) 



Fig. 6. — Satellite location angle (6) for different subsamples of host galaxy fracDeV with 
Icr Poisson error bars. The upper panel provides the following subsample-specific data: 
the number of galaxies (Ngais), the fractional slope as a percent decrease (AA^[0 : 90]), 
uncertainty in the percentage slope (ctaat), and anisotropy significance (IAI/cta). The inset 
directly compares the normalized fractional slopes of the different subsamples. N/Nq is 
normalized galaxy density, where Nq is the galaxy density along the host galaxy projected 
major axis N{6 = 0). The dashed black line indicates zero slope (isotropy). 
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Host Major Axis Preference: satellite galaxy concentration 
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Fig. 7. — Satellite location angle {9) for different subsamples of satellite galaxy fracDeV 
(concentration). The format of this figure is identical to that of FigurelHi AH four subsamples 
show some preference for the host major axis, de Vaucouleurs profile-dominated satellite 
galaxies {fracDeV > 0.5) prefer the major axis most strongly. 
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Host Major Axis Preference: galaxy colors 
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Fig. 8. — Satellite location angle {9) as a function of host galaxy (HG) and satellite galaxy 
(SG) colors. The format of t his figure is identic al to that of Figure [61 We perform color sep- 
aration using the division of iBailin et al.l (120071 ) (Fig. [2]). Major axis preference is restricted 
to satellite galaxies of red hosts and is strongest specifically for red satellites of red hosts. 
Satellites of blue host galaxies are generally consistent with isotropy. 
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Host Major Axis Preference: velocity separation (km/s) 
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Fig. 9. — Satellite location angle (6) for different subsamples of host-satellite velocity sepa- 
ration {AV). The format of this figure is identical to that of Figure [61 Major axis preference 
is restricted to the two lower AV^ bins (AV^ < 600 km/s). At large AV (> 900 km/s), the 
satellite galaxy distribution is consistent with isotropy. 
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Major Axis Preference: radial separation (h"'' l<pc) 
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Fig. 10. — Satellite location angle {9) for different subsamples of projected host-satellite 
radial separation {AR). The format of this figure is identical to that of Figure [61 Major axis 
preference is strongest and most significant at close separations (Ai? < 250 kpc). Despite 
increasing subsample size, this anisotropy decreases steadily in strength and significance with 
increasing radial separation. 
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Fig. 11. — Satellite location angle {6) for different galaxy group sizes (satellites per system, 
Nsg)- The format of this figure is identical to that of Figure El Most satellites exist in 
small systems {Nsg ^ 5). Alignment strength (fractional slope) increases with increasing 
group size. The subsample with the largest groups (21 < Nsg ^ 1000) has the highest trend 
significance despite having the fewest member objects, suggesting that major axis preference 
is more prevalent in larger galaxy groups. 



-45- 



Isophotal Location Angle (6) 



cn 

c 

o 

u 

"5 

c 

o 

-I—' 

u 
o 



E 
o 

cn 
o 

(/) 

X 



SGFD < 0.1 -0.1 < SGFD < 0.5 -0.5 < SGFD < 0.9 - SGFD > 0.9 





I 
o 

V 

o 




p 

A 

o 
■n 

A 
O 




0.5 1 0.5 1 0.5 1 0.5 1 

Radial separation (AR, h~'' Mpc) 

Fig. 12. — Differential distribution of satellite galaxy (SG) angular position (6) relative to 
host galaxies (HGs). Each panel column contains a different range of satellite galaxy fracDeV. 
Similarly, each panel row hosts a different range of host galaxy fracDeV. Projected radial 
separation {AR < l/i~^Mpc), in 10 bins of width 100/i~^kpc, spans the horizontal axis. 
See Fig. [5] for an example of this plotting technique. High significance (fractionally small 
error bars) detections of SG preference for HG major axes (values below 0) are markedly 
concentrated towards the upper-right (high fracDeV). All objects have velocity separation 
AV < 500km/s and good agreement between isophotal and model position angles (APA < 
15°). 
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Hostward Elongation: satellite galaxy concentration 
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Fig. 13. — Hostward elongation (</>) for different subsamples of satellite galaxy fracDeV 
(concentration). The format of this figure is identical to that of Figure El The observed 
hostward elongation trend is restricted to high-/rac-DeF (> 0.9) satellite galaxies. 
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Hostward Elongation: host galaxy concentration 
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Fig. 14. — Hostward elongation (0) for different subsamples of host galaxy /racDe F (concen- 
tration). The format of this figure is identical to that of Figure [61 All four subsamples are 
statistically consistent with each other and with the total group sample. Hostward alignment 
strength thus appears to be independent of host galaxy concentration. 
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Hostward Elongation: galaxy colors 
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Fig. 15. — Hostward elongation (0) as a function of host and satellite galaxy colors. The 
format of this figure is identical t o that of Figure [61 We perform color separation using 
the division of iBailin et al.l (120071 ) (Fig. [2]). Hostward elongation is restricted to the two 
subsamples with red satellite galaxies. Host galaxy color has no discernible effect on hostward 
elongation trend strength. 



- 49 - 



Hostword Elongation: velocity separation (km/s) 
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Fig. 16. — Satellite galaxy hostward elongation (</>) as a function of velocity separation (AV). 
The format of this figure is identical to that of Figure O Each AV subsample is statistically 
consistent with the total group sample. The subsamples are also loosely consistent with 
each other. The high significance in the lowest velocity bin (|Ay| < 300 km/s) is a result 
of subsample size. On its own, host-satellite AV^ does not influence the degree of hostward 
alignment. 
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Hostward Elongation: radial separation (h ^ kpc) 
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Fig. 17. — Satellite galaxy hostward elongation (0) as a function of projected radial sep- 
aration (AR). The format of this figure is identical to that of Figure O No subsamples 
show anisotropy at very high significance. Hostward alignment is strongest in the innermost 
{AR < 250 h-^ kpc) radial bin. 
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Hostward Elongation: nunnber of satellites 
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Fig. 18. — Satellite location angle as a function of the number of satellites per host (i.e., 
group size). The format of this figure is identical to that of Figure [61 All four group size 
subsamples are consistent with other and with the total group population. There is no 
obvious dependence of hostward elongation on group size. 
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Fig. 19. — Holmb erg-like effect seen among select satellite galaxies (SGs). In the above 
histograms we include only groups with inclined {qmom < 0.53 and qiso < 0.53), blue host 
galaxies (HGs) with fewer than 25 identified satellites within 1 Mpc. We further require 
that isophotal and model HG position angle (PA) agree well (APA < 15°). Within these 
systems, we keep SGs that are 15x less luminous than their hosts and satisfy |AV| < 
500 km/s. We eliminate all SGs within 5 Petrosian radii of the host to exclude mis-detected 
bright knots in the host galaxy. The red curve above includes objects with Ai? < 500 kpc 
and Ay < 500 km/s while the blue curve extends only to Ai? < 300 /i"^ kpc and A^ < 
300 km/s. 
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Fig. 20. — Satellite galaxy (SG) distribution {6) trends of the same hosts and satellites 
using isophotal position angle (blue) and galaxy model PA (red). We exclude round objects 
{qiso > 0.9 or qmom > 0.9) for which PA may be inaccurate. The isophotal slope error 
bars, comparable in magnitude to those of the de Vaucouleurs slopes, have been omitted for 
clarity. Comparing these two cases, we observe little difference for any profile type or radius, 
suggesting that the HG population is unaffected by discrepant PAs. 



-54- 



isophotal Hostward Elongation (0) de Vaucouleurs 



0,5 


-0,5 



0,5 


-0,5 



0,5 


-0,5 



0,5 



-0,5 



SGFD < 0.1 



0.1 < SGFD < 0.5 



0.5 < SGFD < 0.9 



SGFD > 0.9 










o 

A 

I 
O 
■n 
o 

A 

O 

to 



o 
O 
A 
O 




ooooo^ooooo^ooooo^ooooo 
NjVoibo ho-i^aico Nj^aico nj^o^cq 

Radial separation (AR, h"'' Mpc) 



Fig. 21. — Hostward elongation (0) in different bins of host galaxy (HG) fracDeV, satellite 
galaxy (SG) fracDeV, and projected radius (AR). The isophotal slope error bars, comparable 
in magnitude to those of the de Vaucouleurs slopes, have been omitted for clarity. We observe 
a marked difference between (j)iso (blue) and (pdev (red) among the same high fracDeVSGs 
(rightmost column). In particular, we observe a striking disparity between isophotal and de 
Vaucouleurs hostward elongation among highly concentrated {fracDeV > 0.9) SGs of highly 
concentrated {fracDeV > 0.9) HGs (top-right panel). We examine this subsample in greater 
detail in Fig. [22j The discovery of this discrepancy motivated our APA < 15° requirement. 
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Fig. 22. — Closer inspection of hostward elongation (0) in the subsample of host-satellite 
pairs with HG fracDeV> 0.9, SG fracDeV> 0.9, \AV\ < 500km/s, and AR < lOO/i-^kpc, 
in which we find the greatest discrepancy between results based on isophote- and de Vau- 
couleurs model-based PAs (see Fig. |2T]). We further separate the galaxies with good PA 
agreement (APA < 15°, blue) from those with poor agreement (APA > 15°, red). The 
whole subsample (no APA requirements, black) is the sum of the red and blue histograms. 
In the figures above, we have excluded all satellites falling within 1 Petrosian radius (rp) 
of the host galaxy which are very likely misidentified bright knots within the host itself. 
Without the APA < 15° requirement, the hostward elongation signal from isophotal PA 
fractional slope is 4 times that we obtain from de Vaucouleurs model PAs. 
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Fig. 23. — The exponential and de Vaucouleurs model-fitting procedures have several known 
flaws. Above, the distribution of position angle (PA) is shown for exponential model flts 
(left) and isophotes (right). While the isophotal PAs have an effectively random distribution, 
the model-based PA suffers from signiflcant anisotropy. For reference, we provide the Ax^ 
between the observed distribution and an isotropic (constant) model. 



Fig. 24. — Spiral structure causes isophotal twisting, which can lead to highly discrepant de 
Vaucouleurs and isophotal position angles. Above, the galaxy models seek out the light in 
the bulge and/or bar, whereas the 25 mag / sq. arcsec isophote roughly detects where the 
spiral arms fade into the sky (on much larger scales). 
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Fig. 25. — Position angles measured with isophotes and galaxy models can provide highly 
discrepant results. Unlike the previous example (Figure [2^ where two discrepant position 
angles both correspond to obvious galactic features, the objects above have nearby neighbors 
which distort the appearance of the isophotes nearby. In several of the cases above, this 
systematic effect could easily mimic a hostward elongation signal. 



